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Abstract 



Three-dimensional spin models of the Ising and XY universality classes are studied 
by a combination of high-temperature expansions and Monte Carlo simulations applied 
to improved Hamiltonians. The critical exponents and the critical equation of state are 
determined to very high precision. 
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1 Introduction 



The notion of universality is central to the modern understanding of critical phenomena. It is 
therefore very important to compare high-precision theoretical and experimental determinations 
of universal quantities, such as critical exponents or universal amplitude ratios, for systems 
belonging to the same universality class. 

There exist several different methods for determining critical quantities. One may study 
lattice models by means of high-temperature (HT) expansions or Monte Carlo (MC) simula- 
tions, or may consider continuum models and apply the well-known methods of perturbative 
field theory. All these methods require some extrapolation: in HT studies one wishes to de- 
termine the behavior at the critical point (3 C from a perturbative series around (3 = 0; in MC 
studies that use finite-size scaling methods, an extrapolation L — > oo is needed; in perturba- 
tive field-theory calculations one must compute the value for g = g*, g* being the fixed point, 
from a perturbative series in powers of g. The accuracy of these extrapolations depends on 
the analytic behavior of the considered functions at (3 C or g*, which, in turn, is determined 
by the renormalization-group (RG) theory. The complex structure of the critical behavior, 
characterized by a multitude of subleading exponents u^, gives rise to nonanalyticities at the 
critical (or fixed) point that make the extrapolations difficult and often introduce large and 
dangerously undetectable systematic errors. We will not review here the field-theoretical case 
and we address the reader to Refs. |l|-f§] and we will only discuss HT and MC methods. 

The precision of the results which can be extracted from the analysis of HT series and 
from MC simulations is mainly limited by the presence of confluent corrections with noninteger 
exponents. Let us consider, e.g., the magnetic susceptibility x(L,T) in a finite box L d , as a 
function of the reduced temperature t. The analyses of HT expansions aim to determine the 
infinite- volume asymptotic behavior for t — > 0, which is 

x{oo,t) = Ct" 7 (l + a 0il t + ...+ a 1A t A + a 1<2 t 2A + ... + a 2il t Aa + ...), (1) 

where A, A 2 , ... are universal exponents. Analogously, MC simulations may try to determine 
the volume dependence of x{L, t) at the critical point (or, in more sophisticated and efficient 
approaches, for a sequence of temperatures t(L) approaching t = as L — > oo): 

X (L,0) = CIP/"(1+ a hl L~ A / u + a lt2 L~ 2A /» + ... + a 2:1 L~ A ^ + ...). (2) 

In both calculations the corrections with exponents A or A/u are one of the most important 
sources of systematic errors. To overcome these problems one may use improved models, that 
is models for which the leading correction to scaling vanishes, i.e. a^i = a^i = 0. Such models 
can be determined by considering a one-parameter family of theories belonging to the given 
universality class, depending, say, on A, and by tuning the irrelevant parameter A to the special 
value A* for which a^i = Si i = 0; we will call such models "improved". 

MC algorithms and finite-size scaling techniques are very effective in the determination 
of A* and /3 C , but not as effective in the computation of critical exponents or other universal 
quantities. On the other hand, the analysis of HT series is very effective in computing universal 
quantities, but not in computing A* and f3 c . 

The strength of the two methods can be combined by computing A* and (3 C by MC, and 
feeding the resulting values into the analysis of HT series (by "biasing" the analysis); this 
greatly improves the quality of the results. 
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Ising 






XY 






7 


V 


V 


7 


V 


a 


IHT* 


1.2372(3) 


0.6301(2) 


0.0364(4) 


1.3177(5) 


0.0380(4) 


-0.0146(8) 


IHT 


1.2371(4) 


0.6300(2) 


0.0364(4) 


1.3179(11) 


0.0381(3) 


-0.0150(17) 


HT 


1.2375(6) 


0.6302(4) 




1.322(3) 


0.039(7) 


-0.022(6) 


MC 


1.2367(11) 


0.6296(7) 


0.0358(9) 


1.3177(10) 


0.0380(5) 


-0.0148(15) 


FT (a) 


1.2396(13) 


0.6304(13) 


0.0335(25) 


1.3169(20) 


0.0354(25) 


-0.011(4) 


FT(b) 


1.2403(8) 


0.6303(8) 


0.0335(6) 


1.3164(8) 


0.0349(8) 


-0.0112(21) 



Table 1: Critical exponents of the three-dimensional Ising and XY models. 



In order to keep systematic errors under control, one may consider several different families 
of models in the same universality class and check that they give compatible results for universal 
quantities. 



2 Critical exponents 

Without further discussion, we present in Table [| a selection of results for the critical exponents 
7, v, and rj of the three-dimensional Ising model; for other exponents, see Ref. ||. IHT denotes 
the results of Ref. ||, where three different improved models were considered. IHT* is a new 
determination in which we bias the analyses by using the MC estimate of /3 C , as we did in 
Ref. || for the XY model. HT is a "traditional" HT determination |7[] obtained by analyzing 
the 25th-order series for the Ising model obtained in Ref. |§ by means of biased approximants; 
MC are Monte Carlo results for the 4 improved Hamiltonians pfl (see also Ref. ||10||); FT 



are results from the field-theoretical expansion in fixed dimension: results (a) are taken from 
Ref. [TT], results (b) from Ref. \T2\. Other results can be found in Ref. ||13|| . The agreement 



among the different determinations is overall satisfactory, although small systematic deviations 
are observed between the lattice (IHT and MC) results and the field-theoretic estimates. We 
suspect that the error estimates of Ref. |L2| are quite too optimistic. 

Similar techniques can be applied to the XY model, with results of comparable quality. We 
present results for the critical exponents 7, 77, and a (we remind that dv = 2 — a) in Table |2j. 
They are taken from Refs. @ (IHT), § (IHT*), p[ (HT), § (MC), (FT(a)), and p 



(FT(b)). Other results can be found in Ref. ||13|| . These results should be compared with the 
precise experimental estimate [16[] a = —0.01056(38) obtain from a Space Shuttle experiment 
for the A transition of 4 He (cf. footnote 2 in Ref. for discussion of the experimental results). 
There is disagreement between IHT* and experiment; it would be interesting to improve further 
the theoretical computation, and to have an independent confirmation of the experimental 
measurement. 
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3 Critical equation of state 



The equation of state relates the thermodynamical quantities M, H, and T. It is easily acces- 
sible experimentally and thus it is important to have predictions for its behavior in the critical 
limit. 

In order to determine the critical equation of state, we start from the effective potential 
(Helmholtz free energy) 

T{M,t) = MH-^\ogZ{H,t). (3) 
In the critical limit J-(M, t) obeys a general scaling law: 

= F{M) - F(0) = t du F sing (MrP). (4) 
In the HT phase, JF can be expanded in powers of M 2 around M — 0. We can write 

AT = JF(M) — JF(0) = —A(z), (5) 

■4W = ^ + ^ + T,^y^, (6) 

where z oc Mt~^ — the normalization is fixed by Eq. (f|) — , m is the second-moment mass, 
and #4 is the renormalized zero-momentum four-point coupling constant. In the critical limit, 
t — > 0, M — > at z fixed, the function A(z) is universal. The (universal) critical limit of and 
r2j can be computed from the HT expansion of the zero-momentum 2j-point Green's functions. 
For the Ising model, we obtain || 

9i = 23.54(4), r 6 = 2.048(5), 
r 8 = 2.28(8), r 10 = -13(4). 

By using (|6|), the equation of state can now be written as 
Equivalently, one can write 

H(M,t) = aM s f(x), (8) 

where x oc tM _1///3 is normalized so that x = — 1 corresponds to the coexistence curve and a is 
fixed by the normalization condition /(0) = 1. The advantage of this representation is that by 
varying x for x > —1 one obtains the full equation of state, while, by using F(z) an analytic 
continuation in the complex plane is needed to reach the coexistence curve. The analyticity 
properties of F(z) and f(x) are constrained by Griffiths' analyticity. 

It is possible to implement all analyticity and scaling properties of the critical equation of 
state introducing a parametric representation [|T7|-|TT| 

M = moR^e, 

t = R{i-e 2 ), 

H = h R p5 h{6), 
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^0 


u 2 


Qc 


Ue 


IHT 


0.530(3) 


4.77(2) 


0.3330(10) 


1.961(7) 


HT+LT 


0.523(9) 


4.95(15) 


0.324(6) 


1.96(1) 


MC 


0.560(10) 


4.75(3) 


0.328(5) 


1.95(2) 


MC 


0.550(12) 








FT 


0.540(11) 


4.72(17) 


0.331(9) 


2.013(28) 



Table 2: Universal amplitude ratios for the three-dimensional Ising model. 



where h(9) is normalized as h(9) = 9 + 0(9 3 ). Note that 9 = corresponds to the HT phase 
t > 0, M = 0, 9 = 1 to the critical isotherm t = 0, and 9 = 9 , where 9q is the first positive 
zero of h(9), to the coexistence curve. The analytic properties of the equation of state are 
reproduced if h(9) is analytic in the interval [0,9 ). Given h{9) the equation of state is easily 
obtained: 

Approximate representations of the equation of state are obtained by approximating h(9) be 
an odd polynomial in 9, i.e. 

h(9) = 9 + Y,h 2 n+i9 2n+1 . (10) 

n=l 

In Ref. the coefficients were determined by using IHT results for the constants r 2n and a 
variational condition. More precisely, one first uses the (k — 1) estimates of r$, r 2 k+2 to 
fix /i 5 , h 2 k+i in terms of h 3 — in Ref. f| the variable p 2 = 6(7 + ^3) was used — and then 
requires physical results to be stationary with respect to variations of h 3 . The idea behind 
this method is that in the exact case the function h(9) is not uniquely defined: There exists 
a one-parameter family of equivalent h{9). Thus, one parameter in h(9) can be fixed at will. 
Whenever we approximate h(9) this is no longer true. What we can require is that physical 
results have the weakest possible dependence on the parameter. 

We use the values of /3, 5, r 6 , r 8 , r 10 obtained by IHT to compute successive approximations 
to h(9); we check the stability of the values of several universal amplitude ratios in order to select 
the best approximation. In Table |3] we report the results of Ref. [||] for the following amplitude 
ratios: U = A + /A~, U 2 = C+/C~, Q c = B 2 {f + f/C + , = /+ / 'f~ . The amplitudes are 
defined in terms of the critical behavior for t — > ± of the specific heat Ch = A \t\~ a , of the 
second-moment correlation length £ = f ± \t\~ u , of the susceptibility \ = C^il -7 , and of the 
spontaneous mag netization M = B(-t) 13 . We compare these estimates with results obtained 
using different methods. HT+LT is a combination of HT and low-temperature expansions 



20|. pl[] ; the other theoretical determinations are the same discussed for the critical exponents, 



and are taken from Refs. (IHT), |22|,[^| (MC), and [24-26[| (FT). The agreement among the 
different determinations is again satisfactory. 

In the XY case, there are Goldstone singularities at the coexistence curve. In three di- 
mensions, the leading singular behavior is correctly reproduced if h{9) ~ {9 — #o) 2 - We may 
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therefore set 

KB) 

The constant 6q and the k coefficients c n are fixed |||I4J by requiring the approximation to 
reproduce the {k + 1) parameters r 6 , . . . , r 2 k+i- 

Only the ratio Uq = A + /A~ is measured experimentally to high precision. [0 Such a 
ratio has been determined in Ref. || using the equation of state and the IHT* estimate of a: 
They obtain Uq = 1.062(4). This result disagrees with the experimental estimate Uo = 1.0442. 
However, the estimate of Uo is strongly correlated to that of a. Indeed, by using a slightly lower 
value of a, a = —0.01285(38), Ref. [pTj] found U = 1.055(3). Thus, the disagreement between 



the estimate of Ref. || and experiment can be reconduced to the discrepancy in a. Other 
estimates of a have been obtained by using the fixed-dimension expansion, [fj4j U = 1.056, 
and the e-expansion, fiZJR U = 1.029(13). 



4 Conclusions 

The study of HT series of "improved" models, with parameters determined by MC simulations, 
allowed us to compute with high precision the universal quantities — critical exponents and 
equation of state — characterizing the critical behavior of the symmetric phase. 

Suitable approximation schemes allow the reconstruction of the critical equation of state 
starting from the symmetric phase; many universal amplitude ratios can be computed. 

For the Ising universality class, theoretical computations are much more precise than ex- 
periments. On the other hand, for the XY class, some very precise experimental results for 
a and A + / A~ have been obtained. [16| There is disagreement with the most precise theoreti- 



cal results. |6| A new-generation experiment is in preparation p8[ ; it would be interesting to 
improve further the theoretical computations as well. 
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